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Abstract 



The effective interaction between a quark and an anti-quark as obtained previ- 
^ ' ously with by the method of iterated resolvents is replaced by the so called up-down- 

Cp . model and applied to flavor off-diagonal mesons including the pion. The only free 

«) , parameters are those of canonical quantum chromo-dynamics (QCD), particularly 

the coupling constant and the masses of the quarks. 

The so obtained light-cone wave function can be used to calculate the pion's 
. form factor, particularly its mean-square radius can be computed analytically. The 

I results allow for the exciting conclusion that the pion is built by highly relativistic 

Q^' constituents, in strong contrast to composite systems like atoms or nuclei with non- 

• relativistic constituents. 

1 Introduction 

• One of the most urgent problems in contemporary physics is to compute the structure of 
^ hadrons in terms of their constituents, based on a covariant theory such as QCD. 

Among the hadrons the pion is the most mysterious particle. I have proposed an 
oversimplified model, the t|-model, which has many drawbacks but the virtue of being 
inspired by QCD and of having the same number of parameters one expects in a full 
theory: namely the 6 flavor quark masses, the strong coupling constant (7) and one 
additional scale parameter (8) originating in the murky depth of renormalization theory. 

The model is QCD-inspired by virtue of the fact, that it is based on the full light-cone 
Hamiltonian as obtained from the QCD-Lagrangian in the light-cone gauge, with zero- 
modes disregarded. In consequence, the pion is treated on the same footing as all other 
pseudo-scalar and pseudo-vector mesons. 

The model should be contrasted to Lattice Gauge Calculations, see for example 
It is not generally known that LGC's have considerable uncertainty to extrapolate their 
results down to such light mesons as a pion. It is also not generally known that lattice 
gauge calculations get always strict and linear confinement even for QED, where we know 
the ionization threshold. The 'breaking of the string', or in a more physical language, the 
ionization threshold is one of the hot topics at the lattice conferences 0. Moreover, in 
order to get the size of the pion, thus the form factor, another generation of computers is 
required, as well as physicists to run them. 
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The model should be contrasted also to phenomenological approaches. They usually 
do not address to get the pion. For the heavy mesons, where they are so successful 0, 
phenomenological model have quite many parameters, in any case more that the above 
canonical ones. A detailed comparison and systematic discussion of the bulky literature 
can however be postponed, until we are ready to solve the full Eg. (pH]) . 

The model should be contrasted, finally, to Nambu-Jona-Lasinio-like models which are 
so successful in accounting for isospin-aspects. I cannot quote the huge body of literature 
but mention in passing that the NJL-models are not renormalizable, that NJL has no 
relation to QCD, and that NJL deals mostly with the very light mesons. There is no way 
to treat the heavy flavors, see also 0]. 

2 Motivation 

The light-cone approach to the bound-state problem in gauge theory p aims at solving 
HlcI"^) — ^^1^)- If one disregards possible zero modes and works in the light-cone 
gauge, the (light-cone) Hamiltonian Hic is a well defined Fock-space operator and given 
in 0. Its eigenvalues are the invariant mass-squares of physical particles associated 
with the eigenstates |\E'). In general, they are superpositions of all possible Fock states 
with its many-particle configurations. For a meson, for example, holds 

l^meson) = T, qq{Xi, k±^, Xi)\qq) + J2 gg{Xi, k±^,, \i)\gg) 
i i 

+ E ^qqg{xi,k^„Xi)\qqg) + E ^qqqq{xi,k±^,Xi)\qqqq) + .... 

i i 

If all wave functions like \E'gg or ^ gg are available, one can analyze hadronic structure in 
terms of quarks and gluons 0. 

For example, one can calculate the space-like form factor of a hadron quite straight- 
forwardly. As illustrated in Fig. |1], it is just a sum of overlap integrals analogous to the 
corresponding non-relativistic formula [§]: 

n<l') = E E / n ^^^^^ ^"1- ^0 ^i^H^. k^^. A.). (1) 

n,Xi o- i 

Here Cq is the charge of the struck quark, A^ ^ g^, and 

^ _ 1 k±i — Xiqj_ + q± for the struck quark 
[ k±i — Xiq_L for all other partons, 

with ql = Q'^ = All of the various form factors of hadrons with spin can be obtained 
by computing the matrix element of the plus current between states of different initial 
and final hadron helicities. 

3 The method of iterated resolvents 

Because of the inherent divergencies in a gauge field theory, the QCD-Hamiltonian in 
3+1 dimensions must be regulated from the outset. One of the few practical ways is 
vertex regularization ^, where every Hamiltonian matrix element, particularly those 
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Figure 1: Calculation of the form factor 
of a bound state from the convolution of 
light-cone Fock amplitudes. The result 
is exact if one sums over all ipn- 



Figure 2: The Hamiltonian matrix for a me- 
son. The matrix elements are represented by 
energy diagrams. Only vertex diagrams V are 
shown. Zero matrices are marked by a dot (■). 



of the vertex interaction (the Dirac interaction proper), is multiplied with a convergence- 
enforcing momentum-dependent function. It can be viewed as a form factor P]. The 
precise form of this function is unimportant here, as long as it is a function of a cut-off 
scale (A). 

Perhaps one can attack the problem of diagonalizing the (light-cone) Hamiltonian 
Hlc by DLCQ, see for example 0. But, alternatively, it might be better to reduce 
the many-body problem behind a field theory to an effective one-body problem. The 
derivation of the effective interaction becomes then the key issue. By definition, an 
effective Hamiltonian acts only in the lowest sector of the theory (here: in the Fock space 
of one quark and one anti-quark). And, again by definition, it has the same eigenvalue 
spectrum as the full problem. I have derived such an effective interaction by the method 
of iterated resolvents , that is by systematically expressing the higher Fock-space wave 
functions as functionals of the lower ones. In doing so the Fock-space is not truncated 
and all Lagrangian symmetries are preserved. The projections of the eigenstates onto the 
higher Fock spaces can be retrieved systematically from the gg-projection, with explicit 
formulas given in 

Let me sketch the method briefly. Details may be found in [§, §. DLCQ with its 
periodic boundary conditions has the advantage that the LC-Hamiltonian is a matrix 
with a flnite number of Fock-space sectors, which we denumerate by n, with 1 < n < N. 
The so called harmonic resolution K = LP~^/ (27r) acts as a natural cut-off of the particle 
number. As shown in Figure 0, i^' = 3 allows for N = 8, and K = 4 for N = 13 
Fock-space sectors, for example. The Hamiltonian matrix is sparse: Most of the matrix 
elements are zero, particularly if one includes only the vertex interaction V. For n sectors, 
the eigenvalue problem in terms of block matrices reads 

n 

Y^{t\H,,iuj)\j){j\^iuj)) = Eiuj) {tl^iuj)), for2 = l,2,...,n. (2) 
I can always invert the quadratic block matrix of the Hamitonian in the last sector to 



define the n-space resolvent G^, that is 



Gniu;) = 1—-. (3) 

Using Gn, I can express the projection of the eigenfunction in the last sector by 

{n\^{u)) = a„M E(^l^^nHlj) (jI^M), (4) 

and substitute it in Eq.(|^). I then get an effective Hamiltonian where the number is 
sectors is diminuished by 1: 

Hn^,{Lu) = Hn{uj) + Hn{uj)Gn{uj)Hn{uj). (5) 

This is a recursion relation, which can repeated until one arrives at the gg-space. The 
fixed point equation E{ijj) = uj determines then all eigenvalues. 

For the block matrix structure as in Figure with its many zero matrices, the reduc- 
tion is particularly easy and transparent. For K = ?> one has the following sequence of 
effective interactions: 

if8 = T8, H-j = T-j + VG^V, He = n + VG7V, H,=T, + VGeV. (6) 

The remaining ones get more complicated, i. e. 

= T4 + VG7V + VG7VGeVG7V, (7) 

Hs = Ts + VG^V + VGqVG5VGqV + VG^V, (8) 

H2 = T2 + VG3V + VG5V, (9) 

Hi = T1 + VG3V + VG3VG2VG3V. (10) 

For K = 4, the effective interactions in Eq.(P) are different, see for example |@], but it 
is quite remarkable, that they are the same for the remainder, particularly Eq.(|IU]). In 
fact, the effective interactions in sectors 1-4 are independent of K: The continuum limit 
K ^ 00 is then trivial, and will be taken in the sequel. 

In the continuum limit, the effective interaction in the gg-space has thus two con- 
tributions: A flavor- conserving piece Ueff-conser = VG3V and a flavor- changing piece 
Ueff-change = VG3VG2VG3V. The latter cannot get active in flavor-off-diagonal mesons. 
Notice that these expressions are an exact result. 



4 The eigenvalue equation in the qq-spaice 



After some approximations 0, the effective one-body equation for flavor off-diagonal 
mesons (mesons with a different flavor for quark and anti-quark), becomes quite simple: 



ml + kl ml + kl 
X 1 — X 



{x, k_i_;\i, A2 



:ii^ 



37r2 



dx'£k[R{x',k[;A) a 
x(l — x)x'{l — x') 

A 



(Ai,A2|5|A;,A'2) (x',^i;A;,A^ 



Here, is the eigenvalue of the invariant-mass squared. The associated eigenfunction 
ijj = is the probability amphtude {x, k±; Ai, X2\4') for finding a quark with momentum 
fraction x, transversal momentum k± and helicity Ai, and correspondingly the anti-quark 
with 1—x, —k± and A2. The mi and Tn2 are (effective) quark masses and a is the (effective) 
coupling constant. The mean Feynman- momentum transfer of the quarks is denoted by 



g2 = g2(^^ j:^. 4) = -1 [(A;i - k[f + (A;2 - k'^f 
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and the spinor factor S = S{x, k±; x', k[_) by 

(Ai, X2\S\X[, A'2) = K^i, Xi)Yuik[, X[)] [v{k',, X',Mk2, X2)] . (13) 

The regulator function R{x', k^; A) restricts the range of integration as function of some 
mass scale A. I happen to choose here a soft cut-off (see below), in contrast to the previous 
sharp cut-off 0. Note that Eq.(p!T|) is a fully relativistic equation. I have derived the same 
effective interaction also with the method of Hamiltonian flow equations, see [|l^ . 

The effective quark masses mi and m2 and the effective coupling constant a depend, in 
general, on A. In the spirit of renormalization theory they are renormalization constants, 
subject to be determined by experiment, and hence-forward will be denoted by mi, m2, 
and a, respectively. In next-to-lowest order of approximation the coupling constant be- 
comes a function of the momentum transfer, a — > a{Q; A), with the explicit expression 
given in 0. 

5 The ||-model and its renormalization 

It might be to early for solving Eq.(|ll|) numerically in full glory like in Ref.0. Rather 
should I try to dismantle the equation of all irrelevant details, and develop a simple model. 

The quarks are at relative rest, when k± = and x = x = mi /{mi + 7712). For very 
small deviations from these equilibrium values the spinor matrix is proportional to the 
unit matrix, with 

(Ai, A2|S|A;a;,) ~ 4mim2 5a„a; (14) 
for details see ||10[. For very large deviations, particularly for k'^ ^ fc^, holds 



Q'^k'^\ and {]i\S\]i)^2k'^\ (15) 



Both extremes are combined in the "||-model" |10 



5" 4mim2 Am-imo , . . . . A^ , , 

^ = ^^ + 2^^^ + 2i?(A,g), with Ri^Q) = j^2- (16) 

It interpolates between two extremes: For small momentum transfer, the '2' generated by 
the hyperfine interaction is unimportant and the dominant Coulomb aspects of the first 
term prevail. For large momentum transfers the Coulomb aspects are unimportant and 
the hyperfine interaction dominates. 

The model over-emphasizes many aspects: It neglects the momentum dependence of 
the Dirac spinors and thus the spin-orbit interaction; it also neglects the momentum 
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Figure 3: Nine contours an(A) are plot- 
ted versus A/A, from bottom to top with 
n = 4, 3, ■ ■ ■ , —3, —4. The contours are ob- 
tained by M2(a,A) = nA^ + M^. The 
thick contour n = refers to the pion with 
= Ml. Mass unit is A = 350 MeV. 
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Figure 4: The pion wave function is 
plotted versus p/(1.338m) in an arbitrary 
normalization. The filled circles indicate 
the numerical results, the open circles the 
analytical function $a(p)- 



dependence of the spin-spin interaction. But the 2 creates havoc: Its Fourier transform 
is a Dirac-delta function with all its consequences in a bound-state equation. 

Here is an interesting point: One is familiar with field theoretic divergences like the 
effective masses and the effective coupling constant. One is used less to "divergences" 
residing in a finite number 2. They must be regulated also, and renormalized. 

In consequence I replace Eq.dlTl) by 



MV(a;, k± 



X 1 — X 



i){x, k± 



a f dx'dPk'^ ( 4mim2 , 2A^ 



Stt^J Jx{l~x)x'{l-x') \ A2 + Q2 



+ TTTT^]^i^'^k[_), (17) 



where k±) = (x, k±; t, i 

For equal quark masses mi = m2 = m, the eigenvalues depend now on three parame- 
ters, the canonical a and m, and the regularization scale A. The dependence can be quite 
strong as seen in Figure |^. There, the lowest mass-squared eigenvalue is plotted versus a 
and A for the fixed quark mass m = 406 MeV. 

Since A is an unphysical parameter, its impact must be removed by renormalization. 
Recently, much progress was made on this question 0, |12] : Adding to R{A, Q) a coun- 
terterm C{K,Q) and requiring that the sum R{K,Q) = R{A,Q) + C{A,Q), and thus 
M2(A;a,m), be independent of A, determines C{A,Q). One remains with R{A,Q) = 
/x^/ (/x^ + Q^)- In line with renormalization theory, one then can go to the limit A — > oo. 
In turn, /i becomes one of the parameters of the theory to be determined by experiment. 
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Table 1: The calculated mass eigenval- 
ues in MeV. Those for singlet-ls states 
are given in the lower, those for singlet- 
2s states in the upper triangle. 
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Table 2: Empirical masses of the flavor- 
off-diagonal physical mesons in MeV. Vec- 
tor mesons are given in the upper, scalar 
mesons in the lower triangle. 



6 Determining the canonical parameters 

The theory has seven canonical parameters which have to be determined by experiment: 
a, II and the 5 flavor masses rrif (if we disregard the top). How can we determine them? 

The problem is not completely trivial. Let me restrict first to the light flavors. With 
rriu = md = m, one has 3 parameters, and in consequence needs 3 experimental data. The 
pion mass M^- = 140 MeV and the rho mass Mp = 768 MeV do not suffice. One needs a 
third datum, the mass of an exited pion, for example. 

Since the mass of the excited pion vr^ is not known with sufficient experimental pre- 
cision, and since the ||-model might be to crude a model to begin with, I choose here 
rriu = rriii = 406 MeV and M^* = Mp = 768 MeV, for no good reason other than con- 
venience. These assumptions are less stringent than they sound, by two reasons. First, 
the rho has a mass less than 2m and should be a true bound state. Second, the Yukawa 
potential in Eq.(|r^) acts like a Dirac-delta function in pairing theory for example: it pulls 
down essentially one state, the pion, but leaves the other states unchanged. 

We thus remain with the two parameters a and /i. Each of the two equations, 
MQ(a,/i) = and Mf{a,fi) = M^* determine a function a{^). Their intersection 



point determines the required solution, which is a = 0.761 and fi = 1.15 GeV |TT[. These 
differ marginally from the previous analysis |]T0|, with fi = 1.33 GeV, for which Figure ^ 
yields a = 0.6904. Once I have the up and down mass, the strange, charm and bot- 
tom quark mass can be determined by reproducing the masses of the K~ and B,^ 
respectively. The parameters in the t|-model can thus be taken as 

a mu = md , , 

0.6904 1.33 GeV 406 MeV 508 MeV 1666 MeV 5054 MeV. ^ ^ 



7 The masses of the physical mesons 

Solving Eq. ( ]T7| ) with the parameters of Eq. (I) generates the mass^-eigenvalues of all flavor 
off-diagonal pseudo-scalar mesons. They are compiled in Table ^ The corresponding wave 
functions are also available, but not shown here. In view of the simplicity of the model, 
the agreement with the empirical values |]13| in Table ^ is remarkable. The mass of the 



first excited states in Table |I| correlates astoundingly well with the experimental mass of 
the pseudo-vector mesons, as given in Table 0. Notice that all numbers in Tables and 
H are rounded for convenience. 

Since the t|-model in Eq.(P^ does not expose confinement one should emphasize 
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Table 3: Flavor-diagonal mass eigenvalues 
in the FM-model for pseudo-scalar mesons 
with the parameter a = (491 MeV)^. 



Table 4: Flavor-diagonal mass eigenvalues 
in the FM-model for pseudo-vector mesons 
with the parameter a = (255 MeV)^. 



that the difference between the physical meson masses in Table |T] amd the sum of the 
bare quark masses is larger than a pion mass. One could call this a kind of practical 
confinement. 

What about the flavor diagonal mesons?- They cannot be a solution to Eqs. (pJ]) or 
([TTD, since the flavor- changing piece of the full effective interaction can generate matrix 
elements between different flavors. Thus far the precise structure of the flavor changing 
part Ueff-change = VG3VG2VG3V has uot bceu analyzed in detail, because it requires a 
considerable effort. 

Rather, the following flavor-mixing model (FM-model [|1^) has been investigated. In 
the FM-model, the full effective Hamiltonian including its flavor mixing is reduced to the 
lowest //-states, i.e. to 

{ff\H,ff\ff) = {i,fj\T + VG^V + VG^VG2VG^V\iPf,j,) = M]^, 6ff, + a. (18) 

Conceptually, it is important that Mjj-, is the eigenvalue of Eq . (p^Tj) . The flavor-mixing 
matrix element a = {ipff\VG3VG2VG3V\ipf,p) depends on the flavors and could be cal- 
culated with a solution of Eqs.(|ll]) or (plT]). In the crude FM-model, however, it is treated 
as a flavor-independent parameter to be flxed by experiment. For 5 flavors one faces thus 
the numerical diagonalization of a 5 x 5 matrix. 

The parameter a for pseudo-scalar mesons was fltted to the mass of the t]', and for 
pseudo-vector mesons to the $, with the results compiled in Tables |]and^ . Three facts, 
however, one gets for free: First, the vr*^ is degenerate in mass with vr^, as well as the p° 
with p^. That they form isospin-triplets is a non-trivial aspect of QCD. Second, both the 
77-77' and the 0;-$ splitting are in the right bull park. Third, that the wave functions of 
the 7r°, 77 and t]' have very much SU(3)-character is even less trivial from the point of 
view of QCD. 



8 The wavefunction of the pion 

For carrying out this programme in practice, I need an efficient tool for solving Eq.(|T7p. 



Such one has been developed recently [|T0| . I outline in short the procedure for the special 



case mi = 7712 = m. I change integration variables from x to by substituting 



^i'^z) = 7T H 1 ^(a^) k±) = . k±). (19) 

^ 2Jm^ + kl + kl Jx{l-x) 
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The variables {kz, fci) are collected in a 3- vector p and Eq.(|T^) becomes 



For the present purpose it suffices to restrict to spherically symmetric s-states (f){p) = (f){p) 
and to apply Gaussian quadratures with 16 points. On an alpha work station it takes 
a couple of micro-seconds to solve this equation for a particular case. The resulting 
numerical wavefunction (j){p) is displayed in Figure ^ and compared with 

^aip) = ^f (^l + , with = 1.338 m. (21) 

Such an analytical form is convenient in many applications. For example, the light-cone 
wavefunction ip{x, k±) can be obtained in closed form by Eq. (p!9|) , i.e. 

^ {2x ~ if + kf 

Kf \ Axil — x) m? 

^{x, k^) = A J (22) 

/x(l -x) (^^ (2x - 1)' + kl 

Ax{l — x) p1 

I can use this to calculate the form factor from Eq.(|l|), and thus the exact root-mean- 
square radius fl^, even in closed form with (r^) = —6 c?F2((5^)/c?<5^|q2=o- 

2^ _ 3 34 + 37s2-41s^ + 15s6 + 36(s)(-8- 16s2 + 21s^- 17s6 + 5s«) 

5(s2 - 1)[4 - 4s2 + 3s4 + 36(s)s4(-2 + s^)] ' ^ ' 

with s = m/pa and the abbreviation h{s) = arctan(Vs2 — l)/\/s'^ — 1. The size of the qq 

1 1 
wavefunction is thus (r^) 2 = 0.33 fm, half as large as the empirical value {r'^)ixp = 0.67 fm. 

This completes the goal: I have a pion with the correct mass, and I have an analytic 
expression for its light-cone wave function. Eq.(^2|) could be used thus as a baseline for 
calculating the higher Fock-space amplitudes, as explained in [||. 

It might well be that the wavefunction in Eq.(^) is consistent with Ashery's experi- 
ment [ITH . 



9 Conclusions 

The proposed pion of the t|-model is rather different from the pions in the literature. I 
have found no evidence that the vacuum condensates are important, but I conclude that 
the pion is describable by a QCD-inspired theory: The very large coupling constant in 
conjunction with a very strong hyperfine interaction makes it a ultra strongly bounded 
system of constituent quarks. More then 80 percent of the constituent quark mass is eaten 
up by binding effects. No other physical system has such a property. 

The effective Bohr momentum of the constituents in the pion turns out as pa = 
1.338 m, see Eg . (^I]) . The mean momentum of the constituents is thus 40 percent larger 
than their mass, which means that they move highly relativistically quite in contrast with 

P 



the constituents of atoms or nuclei. No wonder that potential models thus far have failed 
for the pion. One might mention that lattice gauge calculations use all the computer 
power in this world to generate the potential energy of the quarks and then one uses a 
non-relativistic Schrodinger equation to calculate the bound states. 

All this is to be confronted with the present oversimplified t|-model, which however 
has the virtue to calculate the pion and other physical mesons by a covariant and rela- 
tivistically correct theory. To the best of my knowledge there is no other model which 
can describe all mesons quantitatively from the tt up to the T from a common point of 
view, which here is QCD. 
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